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Preface 


The research investigation described herein was conducted as part of the 
Coastal Inlets Research Program (CIRP) under the work unit “Inlet Laboratory 
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(HQUSACE). Program Monitors for the CIRP at HQUSACE are 
Messrs. Barry W. Holliday, John Bianco, and Charles B. Chesnutt. The Program 
Manager is Mr. E. Clark McNair, Coastal and Hydraulics Laboratory (CHL), 
U.S. Army Engineer Waterways Experiment Station (WES), and CIRP Technical 
Manager is Dr. Nicholas C. Kraus, CHL. 


The mission of the CIRP is to conduct applied research to improve USACE’s 
capability to manage federally maintained inlets which exist on all coasts of the 
United States (including Atlantic, Gulf, Pacific, and the Great Lakes regions). 
Objectives are to (a) make management of channels —-design, maintenance, and 
operation—more effective to reduce the cost of dredging, and (b) preserve the 
adjacent beaches in a systems approach that treats the inlet and beach together. 
To achieve the above objectives, CIRP includes work units on short-wave 
modeling, circulation modeling, channels and adjacent shorelines, inlet scour, 
laboratory investigations, field investigations, and technology transfer. 


The study was conducted by CHL personnel, under the general direction of 
Dr. James R. Houston, Director, and Mr. Charles C. Calhoun, Jr., Assistant 
Director, CHL. Direct guidance was provided by Messrs. C. E. Chatham, 
Chief, Navigation and Harbors Division, and Dennis Markle, Chief, Harbors and 
Entrances Branch, CHL. Experiments were conducted by Messrs. William 
Henderson, Computer Assistant, Cecil Dorrell, and Hugh F. Acuff, Jr., Civil 
Engineering Technicians, under the direction of Mr. William C. Seabergh, 
Principal Investigator, Dr. Jane M. Smith, and Messrs. Gordon S. Harkins and 
Michael J. Briggs, all of CHL. Messrs. Wallace Guy and David Daily of the 
WES Information Technology Laboratory provided instrumentation support. This 
report was prepared by Dr. Smith and Messrs. Seabergh, Harkins, and Briggs. 
Word Processing and formatting were completed by Ms. Myra E. Willis, CHL. 
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Conversion Factors, Non-SlI to 
SI Units of Measurement 


Non-SI units of measurement used in this report can be converted to SI units as 
follows: 
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1 Introduction 


Waves in tidal inlets steepen and break on strong ebb currents. Although wave 
breaking at inlet entrances impacts navigation, sediment transport, and wave 
penetration into the inlet, the wave breaking process in the presence of a current is 
poorly understood. 


Background 


Previous laboratory studies of wave breaking on a current include Hedges, 
Anastasion, and Gabriel (1985); Lai, Long, and Huang (1989); and Suh, Kim, 
and Lee (1994). Hedges, Anastasion, and Gabriel developed a limiting spectral 
shape for waves breaking on a current in deep water and tested it with four 
spectra in a wave-current flume. Suh, Kim, and Lee extended the Hedges, 
Anastasion, and Gabriel formula to finite water depths and tested it with nine 
spectra in a flume. In both studies, little of the data is presented, and results are 
given in the form of limiting spectra. Lai, Long, and Huang performed a detailed 
flume experiment of wave-current interaction kinematics in deep water. They 
observed that linear theory predicted kinematics well, if the Doppler shift is 
included; they confirmed deep-water blocking of waves if the ratio of ebb current 
velocity to wave celerity exceeded 0.25; and they observed a downshifting of the 
peak wave frequency for breaking on a strong current. Ris and Holtuijsen (1996) 
used Lai, Long, and Huang’s deepwater breaking data to evaluate breaking 
criteria and found that the whitecapping formulation of Komen, Hasselmann, and 
Hasselmann (1984) under-estimated dissipation. Supplementing this whitecapping 
with the Battjes and Janssen (1978) breaking algorithm gave significantly better 
agreement with the data. In this study, we will evaluate dissipation formulations 
with a new data set that includes shallow to intermediate relative water depths. 


Objectives 


In this report, wave breaking on a current is examined through physical-model 
measurements in an idealized inlet with a steady ebb current. Wave and current 
measurements are used to evaluate wave dissipation models. The goal of the study 
is to determine a dissipation function for wave breaking on a current that is based 
on integrated wave parameters, is applicable for arbitrary water depths, and is 
robust. 
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Scope 


This introduction described candidate wave dissipation formulations and 
selected previous studies of wave breaking on a current. The experiment 
arrangement (including the laboratory facility, the instrumentation, and the 
experimental series) and experiment procedures and data analysis (including 
sequence of events, calibration, sampling, data analysis methods, data format, and 
example plots) are presented in Chapters 2 and 3. Trends in the measured spectra, 
evaluation and development of dissipation formulations, and application of the 
formulations to calculate wave height decay are presented under “Results” in 
Chapter 4. In Chapter 5, conclusions are summarized. 
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2 Experiment Arrangement 


As part of the Coastal Inlets Research Program, a physical model facility was 
created to address research and field problems of tidal inlets. This tool and 
appurtenances necessary to study inlet problems are discussed in this chapter. 


Idealized Inlet Facility 


An idealized inlet was designed to fit in a 46-m- (150-ft)- wide by 99-m- 
(325-ft)- long concrete basin with 0.6-m- (2-ft)- high walls.’ The approach was to 
design an inlet with simplified bathymetry and fairly steep beach slopes so that 
additional features (such as an ebb shoal) could easily be added. Also, it was 
anticipated that a fine sand would be used as both a tracer and as a fully mobile 
bed, which could be placed over the concrete bottom in a thick veneer. A 1:50 
undistorted scale was assumed to determine reasonable inlet dimensions to model. 
However, other scales can easily be assumed to accommodate the study of 
specific processes due to the simplified bathymetry. 


Figure 1 shows the basin area. The ocean-side parallel contours were deter- 
mined by using an equilibrium profile equation from Dean (1977) 


h=Ax°S (1) 


where h is the stillwater depth, x is distance offshore, and A is determined by the 
sediment grain size.” A value of 0.24 was used, as it represented a relatively 
steep beach. The contoured ocean beach slope extends to the 18.3-cm (0.6-ft) 
mean low water (mlw) depth (or the 9.1-m (30-ft) depth when scaled by 1:50) and 
is then linearly transitioned to the basin floor at a depth of 30.4 cm (1.0 ft) (or 
15.2 m (50 ft)) when scaled by 1:50). The inlet throat region converges to a 
depth of 15.2 cm (or scaled to 1:50, 7.6 m (25 ft)) relative to a mlw datum. The 
minimum width is 244 cm across the inlet between mlw contours (or when scaled 
by 1:50, it represents a width of 122 m (400 ft)). Figure 2 shows the inlet throat 
and entrance channel with parallel jetties which have a spacing of 3.66 m 


'A table of factors for converting non-SI units of measurement to SI units can be found on page vii. 


? For convenience, symbols and abbreviations are listed in the notation (Appendix C). 
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Figure 1. Inlet model research facility 
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(12.0 ft) and extend 5.5 m (18 ft) offshore ( or when scaled to 1:50, spacing of 
183 m (600 ft), extending 275 m (900 ft) offshore. 


Based on Froude's model law (Stevens 1942) and the linear scale of 1:50, the 
model-prototype relations in Table 1 were derived. Dimensions are in terms of 
length (J) and time (#). 


Table 1 

Model-Prototype Scale Relations at 1:50 Undistorted Scale 
Model-Prototype 

ees ne Scale Relation 


tara glo ee uel el = 1:50 


Time (tidal and short wave t= = 1:7.07 
period) 


Other scales may be assumed for the bathymetry, so the different scaling relation- 
ships seen in Table 1 would apply. 


The Idealized Inlet Facility is connected to a large sump (volume of 1.98x10° 
(523,000 gal)) for water exchange so that tides may be produced in the facility’s 
ocean to drive tidal currents into and out of the inlet bay. A constant inflow is 
introduced from the sump into the model ocean while a "rolling" gate either reduces 
or increases flow area over an exit pipe into the sump, which causes ocean rise or 
fall, respectively. The rolling gate is regulated by a controller connected to a 
feedback loop comparing actual to desired water level. The two circular shapes in 
Figure 1 are storage tanks each holding 182,000 2 (48,000 gal) water. They can be 
used to simulate a much larger bay area by storing flood tide water and releasing it 
back to the bay to flow to the ocean during ebb flow. Pumps and control valves 
associated with this procedure are located adjacent to the storage tanks. 


A steady-state flow may also be established for ebbing or flooding currents. The 
piping system is shown in Figure 1. Water is either collected (flood flow) or distri- 
buted (ebb flow) through a system of manifolds in the bay which may be adjusted for 
1, 2, or 3 bay channels or a uniform flow across the bay. Water is either released 
(flood flow) or taken from (ebb flow) the ocean headbay to complete the circulation 
which is energized by the pumps located in the upper left corner of Figure 1. 


Either irregular or monochromatic waves were produced by an 80-ft-long, 
unidirectional plunge-type (vertical motion) wave generator (see Figure 1). 
Unscaled wave periods could be varied from 0.5 to about 3 s and wave heights to 
10 cm (at the generator location and for this particular arrangement of the 
generator). Wave angle could be varied for specific tests by moving the generator 
on its castors. 
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Instrumentation and Calibration 


Wave data were collected with electrical capacitance wave gauges which were 
calibrated daily with a computer-controlled procedure incorporating a least-square 
fit of measurements at 11 steps. This averaging technique, using 21 voltage 
samples per gauge, minimizes the effects of slack in the gear drives and hysteresis 
in the sensors. Typical calibration errors are less than 1 percent of full scale for 
the capacitance wave gauges. Wave signal generation and data acquisition were 
controlled using a DEC MicroVax I computer. Wave data were analyzed with a 
DEC VAX 3600. 


Water velocity data were collected with Sontek 2D Acoustic Doppler 
Velocimeters with a side-looking probe which is oriented to collect in- 
formation on the two velocity components in the horizontal plane. Samples were 
collected at 10 Hz, though the instrument makes 250 pings per second and 
averages for each output sample. Accuracy is +0.5 percent of the measured 
velocity, with resolution of 0.1 mm/s and threshold of 0.1 cm/s. The probe 
samples a 0.25-cm° volume located 5 cm from the sensor heads. 


The sensors were placed as seen in Figure 3. An offshore array of wave 
gauges in front of the wave generator measured the waves outside the region of 
currents. In the inlet region a gauge rack was designed to hold both the wave and 
current meters in a colinear manner, with a 0.61-m (2.0-ft) separation between 
alternating sensors. The rack was then moved to other locations for test reruns of 
the same wave and current condition to complete a data set. 


Experiment Series 


Wave conditions for the study were zero-moment wave height H,,,=3.7 and 
5.5 cm, peak spectral period T,=0.7 and 1.4 s, and incident wave direction 
perpendicular to the jetties and at 20 deg relative to the jetties. All waves were 
generated with a TMA spectral form using a gamma value of 3.30. The current 
velocities were 0, 12, and 24 cm/s as determined fom a reference gauge located 
in the center of the channel seaward of the inlet gorge but between the jetties. 
Magnitude of the current decreased approximately 20 percent seaward of the 
jetties (between current meters 3 and 1). Twelve runs were analyzed for the work 
presented here (two periods, two heights, three current speeds, and normal wave 
incidence). Wave parameters considered in this study are summarized in Table 2 
(lab and prototype scale), where U is the average current magnitude, C is the 
wave celerity, and k is wave number. Table 2 gives the target incident wave 
conditions and maximum current in the inlet throat. The wave and current 
parameters cover a wide range of values, which makes the data useful to evaluate 
the wave dissipation formulations for current-induced wave breaking. Each 
experiment run was repeated three times, first with the wave gauge and current 
meter array centered between the jetties and then offset to the left and to the right 
of the center line by 1 m. Waves and currents at the three positions across the 
channel were similar and were averaged. 
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Table 2. Laboratory Wave Parameters 
3.7,5.5cm | 0.7,1.4s | 0,14, 24 | 0.025- 0.25- 0.4- 0.07- 
cm 0.11 0.63 1.4 0.3 
Prototype | 1.85, 2.75m | 4.9,9.9s | 0.10,1.7 | 0.025- | 0O- 0.25- 0.4- 0.07- 
1:50 m/s 0.11 0.45 | 0.63 1.4 0.3 
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Figure 3. Wave and current meter gauge locations 
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3 Experiment Procedures and 
Data Analysis 


This chapter provides information on data collection and analysis. Some 
selected sample plots of data are shown. 


Sequence of Events and Data Sampling 


Following daily calibration of the wave gauges (Chapter 2), the ebb flow was 
first set in the proper range using a flowmeter installed in the pipeline downstream 
of the pump location. The flow setting then was refined by bringing the current to 
its proper level by using a reference current meter in the inlet channel. With a 
stable flow in the channel, current velocities were collected at all locations for 
about 70 s (700 data samples), then the wave generator was turned on and both 
current and wave data were collected. The wave generator was operated for 
1,020 s. Wave gauge sampling rate was 20 Hz, so 20,400 water elevation data 
points were collected at each gauge and 10,200 velocity data samples (10-Hz 
sampling rate) for each sensor were collected during a run. 


Figure 4 compares target and measured spectra at a wave gauge location in 
front of the wave generator. Figure 5 is a snapshot of waves over a portion of an 
experiment at some of the gauges, first for a waves-only experiment test (upper 
two panels), and then a wave-current experiment (lower 2 panels). 


Data Analysis 


The velocity data were analyzed in the time domain. The pre-wave velocity 
record was averaged over the initial 700 samples. Velocity records after the wave 
generator was in operation were examined individually to determine when a 
steady-state average was reached and a long-as-possible record was averaged, 
usually on the order of 8,000 to 9,000 samples (up to 900 s of data). Figure 6 
shows the currents at stations for a particular run. Initially, there is steady state 
flow without waves, then once wave activity has begun, there is a transitory 
period for the adjustment to a near steady state with oscillations superimposed due 
to wave orbital velocities. 
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Figure 4. Comparison of target and measured spectra for wave generation 
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Figure 5. Snapshot of wave records at wave generator and in entrance channel at 
same time for a test without currents (Experiment 58) and with ebb currents 
(Experiment 70) 
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Figure 6. Time series of currents in channel at meter location 5, run 61 
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Wave data were analyzed in two ways. A down-crossing analysis was per- 
formed on the time series of water elevations as well as spectral analysis using a 
Fast Fourier Transform (FFT). For spectral analysis, the fundamental frequency, 
or af, was 0.00098 Hz (1/1020 s). Forty frequency bands were averaged to 
produce the spectral density estimates. The downcrossing analysis produced the 
calculated parameters shown in Table 3. 


Table 3 
Calculated Wave Parameters from Downcrossing Analysis 


For Rayleigh distribution of wave heights, alpha = 2 and 


WEIBULL ALPHA 

beta = 0.5. Truncation of higher wave heights due to break- 
WEIBULL BETA ing increases alpha. 
NO. OF WAVES Number of waves in record 


H (P=0.5) 
(P = probability) 


T (P=0.5) Median wave period 


The FFT or single channel frequency domain analysis was performed over the 
entire 20,400 data points (at=0.05 sec). In the analysis the mean was removed 
and a cosine square taper applied over 10 percent of the data at the beginning and 
end of the data record. The spectral parameters calculated are listed in Table 4. 
Zero-moment wave heights, peak frequencies, and mean velocities for 12 cases 


Median wave height 


Chapter 3 Experiment Procedures and Data Analysis 


13 


14 


| Table 4 
Calculated Wave Parameters from Single Channel Frequency Domain Analysis! 


Peak frequency, CERC method 
NERS acto Ye cele Peak frequency, single band 
NERD. Neh ya Peak frequency, Delft method 
Peak period, CERC method 
Peak period, single band 

Peak period, Delft method 


Average period, calculated as (EMO/EM2) °° 


are given in Appendix A. The waves were averaged across the left, center, and 
right positions at gauges 2 through 7. 
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4 Results 


1-D Wave-Current Interaction 


The dispersion relationship for waves and currents traveling in the same 
direction or directly opposing each other is (Jonsson 1990, and others): 


w - kU = gk tanhkd (2) 


where w is angular frequency, k is wave number, U is current velocity, g is gravita- 
tional acceleration, and d is water depth. The wave period (T = 27/w) is assumed to 
remain constant as the wave propagates from still water onto the current. For the 
situation where U = 0, Equation 2 reverts to the standard dispersion equation. 

Wave blocking (stopping of waves by an opposing current) occurs for relatively 
strong ebb currents for which there is no solution for Equation 2 (|U|>(g7)/(8 7) in 
deep water and |U|>(gd)* in shallow water). The strong current prevents the wave 
from propagating through the channel, and the wave energy is reflected or dissipated 
by breaking. 


Using & calculated from Equation 2, the wavelength is given by: 


2n 


ji, = 
; (3) 
In shallow water, Equation 2 reduces to: 
je 4 
U + fed . 
and Equation 3 becomes: 
L = (U + yg) T (5) 
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The wave height is determined from the conservation of wave action (Jonsson 1990, 
and others): 


re) EC, = W) 


Ox wo 


r 


mie (6) 


where E is wave energy, C,, is relative group velocity of the waves, x is wave 
propagation direction, and @, is relative angular frequency. Equation 6 assumes no 
dissipation due to breaking or bottom friction. The subscript r represents variables 
measured relative to the current, i.e., variables in a coordinate system moving with 
the current. This one-dimensional formulation was developed under the assumption 
of no refraction or diffraction, which is a reasonable assumption for normally 
indicident waves in the idealized inlet. The wave energy is determined from linear 
wave theory as: 


1 
Eas pgH’ (7) 


where H is wave height and p is water density. The relative angular frequency is 
given by: 


w, = ygk tanhkh (8) 


Equation 8 is similar to Equation 2 for the situation of U=0, but its application is 
different. Equation 8 is used to solve directly for w, with the value of k determined 
from Equation 2. The relative group velocity is given by: 


1 © 2kh 
C = os na ] C2 SS 9 
Bo Ts ( a ©) 


Applying Equation 6 between an offshore Region 1 where the current is negligible 
and a Region 2 in the channel (which may have a different depth and a current) 
gives: 
13, GC . 


@ 


a ei (10) 


1 @, 2 


Solving for the wave height in Region 2 gives: 
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(11) 


where C =a/k is the wave celerity. The second expression on the right-hand side of 
Equation 11 for H, is obtained by substituting @, =w,=@,,+k, U;. If Regions 1 
and 2 are both located in shallow water, Equation 11 reduces to: 


ee eee 
[ ; a (12) 
ed 


Breaking Formulations 


= 


There is little information on breaking criteria for wave breaking on a current. 
Most nearshore breaking criteria neglect current and are based on relative water 
depth, defined as the ratio of wave height to water depth, but existing criteria that 
include wave steepness are good candidates for application on a current, ¢.g., 
Miche’s criterion (1951) given by 


Fax = 9-142 L tanhkd (13) 


where H,),,, is the limiting regular wave height,:L is wavelength, k is wave number, 
and d is water depth. The strength of this relationship is that it reduces to a 
steepness limit in deep water and a depth limit in shallow water, thus incorporating 
both limiting factors in a simple form. Battjes and Janssen (1978) applied the 
Miche criterion with the concept of energy dissipation in a bore (LeMehaute 1962) 
for irregular waves in the following form 


jar 
D = -0.25 OQ, f, mae (14) 
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where D is the wave energy dissipation rate, Q, is the percentage of waves break- 
ing (function of the ratio of root-mean-square wave height to H,,,,), and f,, is the 
mean frequency. 


Another class of breaking relationships with possible application to currents are 
the whitecapping formulations used in spectral wave generation models 
(Hasselmann 1974). Two such relationships are: 


k 
De ag Sian < E(@) (15) 


m 


given by Komen, Hasselmann, and Hasselmann (1984) and Komen et al. (1994), 
where a is a coefficient, S is an integrated relative steepness parameter, w,, is the 
mean angular frequency, k ,, is the mean wave number, E is the energy density 
spectrum, @ is angular frequency; and 


2 295 E( fy pies 
Dit cae US 
ta mo 


given by Resio (1987), where € is a coefficient and fis wave frequency. These 
dissipation relationships were developed for waves in the absence of current, but 
they are applied in this study of breaking on a current to give insight about the 
processes. 


Spectra 


Selected spectra from two muns are shown in Figures 7 and 8 (additional spectra are 
given in Appendix B). The peak periods for the runs were 0.7 and 1.4 s, respectively, the 
incident height for both cases was 5.5 cm, and the current was 24 cm/s. The curves are 
labeled for gauges 0, 1, 3, and 6; with gauge 0 offshore (near the wave generator, average 
of wave gauges 1-7 in Figure 3), gauge 1 approximately 3 m offshore of the jetty heads 
(wave gauge 8 in Figure 3), gauge 3 near the jetty heads (wave gauge 10 in Figure 3), and 
gauge 6 approximately 3 m inshore of the jetty heads. The spectra show interesting 
trends. First, Figure 7 (7, = 0.7 s) shows a significant downshifting of the peak 
frequency from the offshore (gauge 0) toward the inlet throat (gauge 6). The peak period 
increased up to 16 percent for U = 24 cm/s, 9 percent for U= 14 cm/s, and varied by 0-6 
percent for no current. Lai, Long, and Huang (1989) reported a similar trend which they 
attribute to nonlinear side band instabilities. The energy dissipated through breaking was 
extracted at the peak frequency and higher, with the slope of the high-frequency tail of 
the spectra remaining fairly constant. This implies that dissipation is related to energy at 
a given frequency or energy is nonlinearly redistributed to maintain the high-frequency 
slope. The energy in the low-frequency end of the spectra increased, most noticeably in 
the cases with longer peak periods (e.g., Figure 8), with and without current. 
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Figure 7. Spectra for 7, = 0.7 s,H,,. = 5.5 cm, and U = 24 cm/s 
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Figure 8. Spectra for 7, = 1.4 sec, H,,, = 5.5 cm, and U = 24 cm/s 
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Calculating Wave Breaking on a Current 


Dissipation rates 


The motivation for these laboratory experiments was to find a wave dissipation 
formulation for wave breaking on a current that can be applied in numerical 
wave-transformation models. Two likely parameters to correlate with dissipation 
are wave steepness (as in deepwater breaking and whitecapping relationships) and 
wave height (as in bore dissipation models). Also, existing dissipation formula- 
tions discussed in the introduction are evaluated with dissipation calculated from 
the measurements. 


Dissipation was calculated from the laboratory measurements by applying the 
action balance equation, Equation 6, modified to include energy dissipation: 


(CG & 
pos is om 


Ox 


The action balance equation was applied between two wave gauges to solve for the 
dissipation D. Figures 9 and 10 show the calculated dissipation as a function of 
steepness and wave height, respectively. Figure 9 shows that the steepness 
parameterization segregates the data by peak wave period. Although dissipation 
increases with steepness, for a given steepness, dissipation is higher for longer 
peak periods. This result foreshadows that the whitecapping dissipation formula- 
tions, which are strongly a function of wave steepness, will not provide good 
estimates of dissipation for this data set. The calculated dissipation is highly 
correlated with wave height, and wave period does not seem to be a controlling 
parameter (Figure 10). This wave-height dependence implies that bore-type 
dissipation formulations, which are functions of wave height, are good candidates 
for estimating dissipation for this data set. Figure 11 shows the Miche limit 
(Equation 13) in terms of maximum wave steepness as a function of relative depth 
(solid line). The measurements, some of which correspond to breaking and some 
do not, also are plotted. The Miche criterion serves as a conservative upper limit 
to the data. The conservatism is not surprising because the formulation is for 
regular waves, and the data correspond to irregular waves. 


Three dissipation formulations were evaluated with the laboratory data, those of 
Komen et al. (Equation 15), Resio (Equation 16), and Battjes and Janssen (Equation 
14). Although current does not enter explicitly in any of these formulations, current 
has been included in the calculation of wavelength and wave number using conser- 
vation of waves and linear theory (see, e.g., Jonsson (1990)). The Resio formula- 
tion is given in Figure 12 as the solid line. This formulation slightly overpredicts 
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Figure 9. Dissipation versus steepness 
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Figure 10. Dissipation versus wave height 
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H/L = 0.142 tanh kd 


Figure 11. Miche criterion 
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Figure 12. Resio dissipation 
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dissipation for the 0.7-s waves, but significantly underpredicts dissipation for the 
1.4-s waves. The correlation coefficient (r*) for the Resio formulation (both 
periods) is 0.46. Likewise, the Komen et al. formulation (not shown) significantly 
underpredicts the dissipation for all wave conditions (7? = 0). Because these two 
formulations are strongly dependent on wave steepness, it is expected that they 
correlate poorly with the data set (see Figure 9). It is not surprising that these 
white-capping expressions do not represent wave breaking on a current in shallow 
water, because the formulations were developed to balance excess wind input in 
saturated spectra without a current. The formulation of Battjes and Janssen is 
plotted in Figure 13 with the laboratory measurements and gives the correct trend in 
the data, but produces wide scatter (large deviation from the formulation), with a 
correlation coefficient of 0.25. 


The poor performance of the dissipation formulations led to development of an 
alternative relationship. Figure 10 showed good correlation of dissipation to wave 
height, so a form consistent with dissipation in a bore was assumed, and the 
following relationship was developed using linear regression: 


1 
D=-0.002 (| 2 (ele lel) store eh Te, (18) 


where H, is a critical wave height below which no dissipation occurs, 


H, = 0.08 L tanh(ka) (19) 


The wavelength and wave number in Equation 19 include modification by a current. 
Equation 18 was chosen to be of the form of Miche’s criterion. Equation 18 is 
plotted against the laboratory data in Figure 14 and shows reasonable agreement 
(r? = 0.78). 


Wave heights 


Applying the action balance equation (Equation 17) between the wave gauges, 
together with the various dissipation models, gives a simple one-dimensional 
shoaling and decay model. The Miche criterion also was applied as a dissipation 
function by limiting the wave energy based on the maximum wave height given by 
Equation 13. This transformation technique was used to calculate wave shoaling 
and breaking through the gauge array for each of the 12 runs. Cases without current 
had little or no wave height decay through the gauge array. Example cases with the 
greatest dissipation are presented in Figures 15-17 (additional cases are given in 
Appendix B). Note in the figures that the second point (x = 120 cm) is consistently 
lower than the nearest points, probably due to a gauge calibration problem. This 
gauge also accounts for some of the scatter in the dissipation calculations. 
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Figure 13. Battjes-Janssen dissipation 
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Chapter 4 Results 


The results generally cluster into two categories. The first category includes the 
formulations of Komen, Hasselmann, and Hasselmann; Miche; and Resio (for 7, = 
1.4 s), which show little or no dissipation for almost all runs and thus significantly 
overpredict the wave height through the gauge array. In Figures 15 and 17, the 
Komen, Hasselmann, and Hasselmann and Miche curves generally overlay each 
other and predict no dissipation. Even in the cases with the highest wave height 
and strongest current (Figures 16 and 18), the Komen et al. formulation predicts 
little or no dissipation, whereas the Miche criterion predicts some wave height 
decay. 


The second category includes the formulations of Battjes-Janssen, the present 
study (Equations 18 and 19), and Resio (for 7, = 0.7 s). These formulations fall on 
the lower portions of the figures and agree well with the measured wave heights for 
the twelve cases, with the exception of cases with T, = 0.7 s and U = 14 cm/s 
(Figure 17). The reason for the poor agreement with this case is unclear. The 
correlation coefficient for the Resio formulation (both peak periods) was 0.71, for 
the Battjes and Janssen formulation 0.86, and for the present study formulation 
0.87. Therefore, in spite of the large scatter between the Battjes-Janssen dissipation 
relationship and the data, the prediction of wave height is generally good, essentially 
equivalent to the expression developed in this study. 
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Figure 15. Wave height transformation for 7, = 1.4. s, H,, = 5.5 cm, and 
U= 14 cm/s 
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Figure 16. Wave height transformation for 7, = 1.4 s, H,, = 5.5 cm, and 
U = 24 cm/s 
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Figure 17. Wave height transformation for 7, = 0.7 s, H,,. = 5.5 cm, and 
U = 14 cm/s 
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Figure 18. Wave height transformation for 7, = 0.7 s, H,, = 5.5 cm, and 
U = 24 cm/s 
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5 Conclusions 


Results from a laboratory study of wave breaking on an ebb current at an 
idealized inlet are presented. The measurements show increased wave shoaling 
and breaking in the presence of the ebb current, as is expected from the conserva- 
tion of waves and linear wave theory. The measurements also show a down- 
shifting of the peak frequency (peak period became longer as much as 16 per- 
cent). Energy was dissipated at the peak and higher frequencies of the spectra. 
Dissipation rates calculated from the measurements were proportional to wave 


height. 


Dissipation algorithms were examined with the data. Whitecapping formula- 
tions (Komen, Hasselmann, and Hasselmann 1984; Komen et al. 1994; Resio 
1987), which are strongly dependent on wave steepness, generally underpredicted 
dissipation. Application of Resio’s whitecapping formulation gave a correlation 
coefficient of 0.46 for predicting dissipation through the idealized inlet. The 
Battjes and Janssen breaking algorithm worked well for predicting wave height 
through the idealized inlet, despite considerable scatter in the dissipation predic- 
tion (7 = 0.25). A relationship for dissipation as a function of wave height 
squared gave improved agreement between calculated and predicted dissipations 
(7 = 0.78), but no substantial improvement over the Battjes and Janssen formula- 
tion for modeling the wave height. 


A dissipation function or breaking criterion applied at a coastal inlet must 
include relative depth and wave steepness, as well as wave-current interaction. 
Depth-limited breaking is more important for longer period waves and steepness 
for shorter period waves. 


The next phase of experiments in the idealized inlet includes construction of an 
ebb shoal and extension of the measurement array farther offshore. Breaking 
over the ebb shoal, wave blocking, scaling relationships, and oblique wave 
incidence will be studied. 
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Photo 5. Wave height: 5.5 cm, wave period 0.7 sec, ebb current: 12 cm/sec 
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Table A2 
Still-Water Depths at Gau 
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Note: Left, center, and right are gauge lines looking seaward from inlet. 
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Notation 


Constant in Equation 1 

Constant in Equation 9 

Wave celerity, m/s 

Wave group celerity, m/s 

Wave group celerity relative to the current, m/s 
Energy dissipation, m7/s 

Water depth, m 

Wave energy or spectral energy density, m?/Hz 
Wave frequency, s 

Mean wave frequency, s 

Gravitational acceleration, m/s” 

Still-water depth 

Wave height 

Critical wave height, above which breaking occurs, m 
Maximum wave height, m 

Zero-moment wave height, m 

Wave number, m” 

Mean wave number, m™ 

Length 

Airy wavelength, cm 

Area, m” 

Volume, m? 

Linear scale of the model 

Percentage of waves breaking 

Steepness parameter 

Time, s 

Wave period of peak energy density of spectrum, s 
Current velocity, m/s 

Volume, m° 

Velocity, m/s 

Coordinate axis, m 

Constant in Equation 16 

Constant 

Water density, kg/m° 
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C1 


C2 


Angular frequency, s” 
Angular frequency relative to the current, s” 
Mean angular frequency 
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